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Abstract: We show that color-kinematics duality is present in tree-level amplitudes of 
quantum chromodynamics with massive flavored quarks. Starting with the color struc¬ 
ture of QCD, we work out a new color decomposition for n-point tree amplitudes in a 
reduced basis of primitive amplitudes. These primitives, with k quark-antiquark pairs and 
(n — 2k) gluons, are taken in the (n — 2)\/k\ Melia basis, and are independent under the 
color-algebra Kleiss-Kuijf relations. This generalizes the color decomposition of Del Duca, 
Dixon, and Maltoni to an arbitrary number of quarks. The color coefficients in the new 
decomposition are given by compact expressions valid for arbitrary gauge group and repre¬ 
sentation. Considering the kinematic structure, we show through explicit calculations that 
color-kinematics duality holds for amplitudes with general configurations of gluons and 
massive quarks. The new (massive) amplitude relations that follow from the duality can 
be mapped to a well-defined subset of the familiar BCJ relations for gluons. They restrict 
the amplitude basis further down to (n — 3)! (2 k — 2)/k\ primitives, for two or more quark 
lines. We give a decomposition of the full amplitude in that basis. The presented results 
provide strong evidence that QCD obeys the color-kinematics duality, at least at tree level. 
The results are also applicable to supersymmetric and D-dimensional extensions of QCD. 
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1 Introduction 


The sheer abundance of QCD background processes at the Large Hadron Collider has 
stimulated remarkable theoretical progress in the treatment of perturbative gauge theory. 
Calculations of previously intractable processes involving high numbers of loops and legs 
have been successfully carried out, such as the inclusive W + 5-jet production at next-to- 
leading order [1], or the inclusive Higgs production at (next-to) 3 -leading order [2]. Mod¬ 
ern methods have been essential for many aspects of this progress, such as the unitarity 
method [3-5], on-shell recursion [6, 7], more recent unitarity-based one-loop methods [8- 
19], and advanced integration techniques [20-26], to mention a few. Even so, increasingly 
refined methods are needed to curb the factorial growth of complexity characteristic of 
perturbative computations. 

The increase of both the number of loops and legs is an unavoidable consequence 
of precision QCD. In the physical observables, contributions from different numbers of 
loops and legs are tied together through the singularity structure and unitarity. High- 
multiplicity tree amplitudes enter explicitly in the real emission contributions, as well as 
secretly in the virtual loop contributions through the modern unitarity-based methods. 
Therefore, progress in treatment of tree amplitudes can be directly or indirectly translated 
to progress in precision QCD phenomenology. 

Modern methods usually decompose tree amplitudes into purely-kinematic primitive 
amplitudes and their coefficients that depend only on color. Such primitive amplitudes 
are gauge-invariant and are often computed via off-shell [27] or on-shell recursion [6, 7], 
which admits an analytic all-multiplicity solution [28, 29] for massless QCD with up to four 
quark-antiquark pairs. The primitives can be assembled in several superficially different 
ways, with varying degree of efficiency, to obtain the full amplitude. The standard SU(IV C ) 
trace-based decomposition of an n-gluon amplitude involves an overcomplete summation 
over (n — 1)!/2 linearly-dependent primitives. The linear dependence comes from two 
sources: from the color algebra of the gauge theory [30] and from an observed kinematic 
algebra, which is closely tied to the former through the color-kinematics duality [31, 32]. 
The corresponding relations satisfied by the primitives are known to all multiplicity: the 
Kleiss-Kuijf (KK) relations [30], and the Bern-Carrasco-Johansson (BCJ) relations [31]. 

More efficient ways of assembling the primitives are known. The Del Duca-Dixon- 
Maltoni (DDM) [33, 34] color decomposition removes all the pure-gluon primitives that are 
redundant under the KK relations, and thus involves only the (n — 2)! members of the KK 
primitive basis. For amplitudes with a single quark-antiquark pair a very similar decom¬ 
position can be used [35, 36]. For general QCD amplitudes, with k quark-antiquark pairs 
and (n — 2k) gluons, much less is known about non-redundant decompositions involving 
well-defined and simple primitive amplitudes. Recently, Melia proposed [37, 38] a basis of 
primitives, for general n and k, that enjoys the same useful properties as the KK basis. For 
distinctly-flavored quarks the basis consists of (n — 2)\/k\ planar color-ordered amplitudes. 
In this paper we give a color decomposition that uses precisely these primitives. Thus we 
obtain a generalization of the DDM decomposition valid for all tree amplitudes in QCD, 
applicable to any gauge group and any representation. 
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The color-kinematics duality [31, 39] highlights the fact that gauge theories generically 
possess an underlying kinematic structure that controls the theory, similar to the way the 
color Lie algebra defines the theory. This structure has been used successfully for tree and 
loop amplitudes in massless gauge theories with and without supersymmetry, including 
massless QCD with no quarks [32, 40-51]. The duality imposes the BCJ relations on pure- 
gluon tree amplitudes that further reduce the basis of independent primitives to one of 
size (n — 3)!. These relations have been observed to still hold in certain specific excursions 
away from the massless pure-gluon case, e.g. for one massless quark-antiquark pair [52, 53], 
or for two or three massive particles [54, 55]. The color-kinematics duality was extended 
to theories with fundamental massless matter in refs. [49, 50], but the corresponding BCJ 
relations were not worked out there. They are instead presented here for matter with any 
mass. 

In this paper, we show that the color-kinematics duality can be straightforwardly 
extended to all tree-level amplitudes in QCD. Furthermore, for amplitudes with k massive 
distinctly-flavored quark-antiquark pairs and (n — 2k) gluons, we find all BCJ relations 
generated by the duality. For k > 2 they reduce the number of independent primitives 
to (n — 3)!(2 k — 2 )/k\, and for k = 0,1 the counting is the standard (n — 3)!. Using this 
reduced BCJ basis, we construct a second new decomposition that uses these primitives 
only. We expect that the BCJ relations obtained here are applicable beyond the case of 
QCD. In particular, super-QCD and D-dimensional QCD should obey the color-kinematics 
duality at tree level just as well. 

This paper is organized in an example-driven way. After we review the color and 
kinematic structure of gauge theory amplitudes in section 2, we proceed to the lower- 
point examples that lead to the main result of section 3 - the new color decomposition 
for a generic QCD amplitude (section 3.3). We return to the lower-point examples in 
the beginning of section 4, where we show that they respect the color-kinematics duality 
and thus satisfy the massive BCJ relations. Then we generalize them and derive a new 
amplitude basis in section 4.5. Our final amplitude decomposition is given in section 4.6. 
We conclude by discussing our results in section 5. 

2 Review and preliminaries 

In this section we review some general properties of the color and kinematic structure of 
tree-level scattering amplitudes in QCD. 1 Amplitudes involving only gluons or at most 
one quark-antiquark pair have a similar form and structure, and are well studied in the 
literature. Adding more quarks makes the composition of the amplitude more involved, 
and this is the topic of the bulk of this paper. 

QCD is a renormalizable gauge theory with only cubic and quartic interactions. In 
what follows we center the discussion around the cubic interactions. The role of the quar¬ 
tic ones is to make the amplitudes constructed from the Feynman rules gauge invariant, 

1 In this paper by QCD we mean Yang-Mills theory with gauge group G and with Nf massive Dirac 
fermions in the fundamental representation (quarks). Additionally, we allow for supersymmetric extensions. 
The cases of Af = 0, 1, 2 (super-)QCD are all included in the general treatment that follows. 
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Table 1: Number of cubic graphs, z/(n, k), in the full n-point amplitude with k distinguishable 
quark-antiquark pairs and (n — 2k) gluons. 


and they carry no new physical information with respect to the cubic interactions. This 
nontrivial statement is made apparent by the on-shell recursion [6, 7], which relies only on 
input from the three-point amplitudes of the theory. 

For the color structures the redundancy of quartic interactions is clear from inspecting 
the four-gluon Feynman vertex 
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( 2 . 1 ) 


where we, for later convenience, use imaginary structure constants f abc = i\/2f ahc . This 
vertex contains the same color factors as the s-, t- and '(/-channel diagrams that are con¬ 
structed from the three-gluon vertices. Hence the quartic vertex can always be absorbed in 
some way into cubic trees without changing the general color structure of a QCD amplitude. 

Therefore, without loss of generality we can write a QCD tree amplitude in terms of 
cubic graphs only. This gives us an expansion of the n-point tree amplitude of the form 2 
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( 2 . 2 ) 


where c,; are color factors, rii are kinematic numerators, and D* are denominators encoding 
the propagator structure of the cubic graphs. The denominators may contain masses, 
corresponding to massive quark propagators. 

Amplitudes with multiple quarks of the same flavor and mass can be obtained from 
distinct-flavor amplitudes by setting masses to be equal and summing over permutations of 
quarks. Therefore, we do not lose generality by taking all k quark-antiquark pairs to have 
distinct flavor and mass. For explicitness, in table 1 we provide total counts of cubic graphs 
for different amplitudes up to eight particles and four quark pairs. This count appears in 
eq. (2.2) as u(n, k). It agrees with the usual counting of standard QCD Feynman diagrams 
restricted to those diagrams that only have trivalent vertices. 

2 Following ref. [31] we absorb all factors of i into the numerators, which is convenient for tree amplitudes. 
The numerators in ref. [50] have a factor of —i pulled out relative to this convention. 
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Figure 1: Color vertices with planar ordering consistent with the color-stripped Feynman rules. 


The kinematic numerators Hi are in general not gauge-invariant and thus have no 
unique expressions. However, certain linear combinations of rii/Di correspond to gauge- 
invariant primitive amplitudes that can be constructed from the color-stripped Feynman 
rules [56] summarized in appendix A. 

The color factors q in eq. (2.2) are constructed from the cubic graphs using only two 
building blocks: the structure constants f abc for three-gluon vertices and generators T°j for 
quark-gluon vertices, as shown in figure 1. When separating the color from kinematics, 
the diagrammatic crossing symmetry only holds up to signs dependent on the permutation 
of legs. These signs are apparent in the total antisymmetry of f abc . Since the structure 
constants can be thought of as the generators in the adjoint representation, (T“ d j)ft c = f bac , 
it is convenient to introduce a similar antisymmetry for the fundamental generators, 

7T = -7T 4^ f cab = _ fbac . ( 2 . 3 ) 

In formulas with suppressed indices we denote the flipped generator by T a = T£j. 

Color factors obey simple relations arising from the Jacobi and commutation identities, 

Jdac Jcbe Jdbc Jcae 

rnCL r T 1 b_ _ rrib /TPO_ 

1 ij 1 jk J 1 jk 

depicted diagrammatically in figure 2. They both imply color-algebraic relations of the 
schematic form 

C/ Cj — Cfc (2.5) 

for the triplets of diagrams ( i , j, k ) that differ only by the subdiagrams drawn in figure 2, but 
otherwise have common graph structure. The interdependence among the color factors c* 
means that the corresponding kinematic coefficients rii/Di are in general not unique, as 
should be expected from the underlying gauge dependence of the numerators. In section 2.2 
we return to this kinematic numerator freedom. Before that, we review how to assemble 
the diagrams into gauge-invariant building blocks, i.e. the primitive amplitudes. 

2.1 Color decomposition of tree amplitudes 

A classic way to remove the relations among color factors is to replace all structure constants 
by generators, 

fabc _ rj-iarpbrjnC^ _ ^rpbrj-iarpC'j ^ gj 

and then eliminate all contracted adjoint indices using the SU(AT C ) Fierz identity, 

T Pii = iah 3 -^-ii 3 6 a - (2.7) 


= f abc f dce , (2.4a) 

= f abc Tl, (2.4b) 
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Figure 2: Color-algebra relations in the adjoint (a) and fundamental representation (b). 
The color-kinematics duality requires that the kinematic numerators satisfy the corresponding 
kinematic-algebra relations, which can be represented by the same graphs. 




This leads to a basis of linearly-independent color structures. In the pure-gluon case it 
gives the familiar color-trace decomposition [36, 57-59], 

A^= Tr(T ai TM 2 ) ...T a -M)A(l,a(2),... ,a(n)), (2.8) 

o-eSn—i({2,...,n}) 


where the sum is over (n— 1)! primitives because the cyclic symmetry of the trace allows one 
to fix the first argument. 3 In a similar but much more cumbersome way, the color content 
of amplitudes with k quark-antiquark pairs and {n — 2k) gluons can be reduced [36, 60-67] 
to color structures of the following type: 


Nl 


(j' a 2k +1 T a h) f'I' a l 1 + l r J' ai A (j lCll k-l + 1 

V ’ /ZlCh V ’ / ioOLO \ * / 


' ^k^k 


(2.9) 


where If € {2k + 1,... ,n}, d/ € {ji,... ,%}■ and (jp + 1) counts the number of disjoint 
cycles in the permutation a. 

An obvious feature of the above SU(AT C ) decomposition is that it is specific to the 
gauge group. A more interesting drawback is that it usually maps the color space to 
a basis that is larger than the number of linearly independent color factors q. Hence 
the resulting kinematic coefficients - the color-ordered amplitudes - are not the minimal 
set of primitives. In other words, they are not independent. Indeed, the purely gluonic 
color-ordered amplitudes can be reduced using the Kleiss-Kuijf (KK) relations [30]. By 
construction, the number of the primitives independent under such color-algebra relations 
must coincide with the number of linearly independent color factors. The KK amplitude 
relations can be written as 


A(l,0,2,a) = (-1)M Y, ^( 1 , 2 ,< t ), ( 2 . 10 ) 

a^aLL\/3 T 

where the sum runs over the shuffle product of the ordered sets a and (3 T , the latter being 
(3 in reverse order. This gives all partially ordered permutations that respect the element 
order of the two sets. The KK relations let us fix the second argument in the pure gluon 
primitives, and hence they reduce the basis to (n — 2)! elements - the KK basis. 

3 The reversal symmetry A( 1, 2,..., n) = (—1 ) n A(n ,..., 2,1) further reduces that count to (n — l)!/2. 
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Figure 3: Multi-peripheral cubic diagram for the color factors in formulas (2.11) and (2.12). All 
permuted legs are gluons, while the horizontal line can be either a quark or a gluon line. 

An interesting exception to the redundancy of the aforementioned decomposition al¬ 
gorithm is the case of the amplitude with a single quark line [35, 36], 

<? = E (™ 3) ■ ■ • ^^) j2n A(l, 2, <r(3),... , a(n )), (2.11) 

a£S n - 2 ({3,...,n}) 

where we have denoted the quark-antiquark pair by a bar below and above the labels, i.e. 
as 1 and 2. The sum is over the basis of (n — 2)! primitives, which are already independent 
under color-algebra relations. Moreover, the color factors in eq. (2.11) are valid for any 
gauge group, as can be guessed from the absence of explicit factors of N c . 

The reason for the nice properties of eq. (2.11) is that this example happens to coincide 
with another basis of color factors that have these properties more generally. For k < 2 the 
color factors in this basis correspond to the “multi-peripheral” graphs shown in figure 3. 
The above k = 1 case can be easily mapped to the pure-gluon configuration by replacing 
T a —> T“ d j. This gives the decomposition of Del Duca, Dixon and Maltoni (DDM) [33, 34]: 

M^ e 0 e = E f a2a ^ bl f bia ^ b2 ...f bn ~ 3a ^ ai A(l,2,a{3),...,a(n)). (2.12) 

cr£S n -2({3,...,n}) 

Note that it is a substantial improvement over the trace decomposition (2.8), since it avoids 
using (n — 2) 2 (n — 3)! primitives altogether. 

The details become more involved when considering generalizations along the lines of 
the DDM decomposition to more than one quark line. A basis of amplitude primitives for 
generic tree amplitudes in QCD, which have the same count as the number of indepen¬ 
dent color structures, was recently found by Melia in refs. [37, 38]. However, it remained 
unknown what are the corresponding color coefficients in a decomposition using that basis. 

In section 3 we discuss the Melia basis in detail, and give the complete amplitude 
color decomposition in terms of this basis. The new decomposition can be thought of as a 
natural generalization of the DDM decomposition to the case of k quark-antiquark pairs, 
in analogy to how the Melia basis is the multi-quark generalization of the KK basis for 
amplitude primitives. Similarly to the multi-peripheral formulas (2.11) and (2.12), the color 
coefficients of the new decomposition will be constructed from the cubic color factors c* 
and will thus be valid for any gauge group (and any group representation for the quarks). 

2.2 Color-kinematics duality 

Let us return to the trivalent graph expansion (2.2) of a gauge-theory amplitude, involving 
kinematic numerator factors n,, color factors c, and denominators D*. As already explained, 
the numerators are not uniquely defined, since, for instance, the quartic-vertex contact 
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terms can be freely absorbed into the cubic graph numerators in more than one way. More 
generally, the ambiguity of the rij is a necessary consequence of the gauge dependence of 
individual Feynman diagrams, and hence the operation of shifting the numerators while 
leaving the amplitude invariant is called generalized gauge transformation [31, 32]. 

It was observed by Bern, Carrasco and one of the current authors (BCJ) [31], that 
within the freedom of shifting the numerators there exists particularly nice choices, such 
that the resulting kinematic numerators rii obey the same general algebraic identities as the 
color factors c*. That is, there is a numerator relation for every color Jacobi/commutation 
relation (2.4) and a numerator sign flip for every color factor sign flip (2.3): 

rii — rij = nk Cj — Cj = Ck , (2.13a) 

rii —» —rii c. L —> —Q . (2.13b) 


Amplitudes that satisfy these relations are said to exhibit color-kinematics duality. While 
both eqs. (2.13a) and (2.13b) should be imposed, the latter is often omitted as the numer¬ 
ators inherit this antisymmetry from the cubic vertices of color-ordered Feynman rules. 

The relations in eq. (2.13) define a kinematic algebra of numerators, which is suggestive 
of an underlying kinematic Lie algebra. While not much is known about this Lie algebra, 
which should be infinite-dimensional due to the continuous nature of momentum, in the 
special case of the self-dual sector of QCD the kinematic algebra is known to correspond 
to certain area-preserving diffeomorphisms [68]. 

A very useful aspect of the color-kinematics duality is that once the numerators satisfy 
eq. (2.13), they can take the place of the color factors in eq. (2.2), since they satisfy the 
same algebraic identities. Through this procedure one obtains a double copy amplitude, 
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tree 

n,k 



n —2 


u(n,k) 
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cubic graphs 


riiHi 

Di 


(2.14) 


which describes scattering in a gravitational theory. More precisely, the scattering ampli¬ 
tude (2.14) will involve (n — 2k) gravitons and 2k matter particles [50]. The tilde notation 
is introduced since the two copies of numerators may be not identical, they can differ by a 
generalized gauge transformation, or the states may differ on the two sides. More generally, 
the two sets of numerators in the double copy construction do not have to belong to the 
same gauge theory [31, 32], producing a wide range of gravity theories with and without 
supersymmetry. For example, pure Yang-Mills theory “squares” to gravity coupled to a 
dilaton and anti-symmetric tensor. Pure Einstein gravity can be obtained by removing 
these extra particles via a ghost-like double-copy prescription for massless quarks [50]. 
In contrast, a highly asymmetric double copy is needed for the amplitudes that couple 
Yang-Mills theory to gravity [69]. 

At tree level the double copy construction is known [31, 39] to be equivalent to the 
field-theory limit of the Kawai-Lewellen-Tye (KLT) relations [70] between open- and closed- 
string amplitudes. However, the color-kinematics duality (2.13) and the double copy (2.14) 
can be argued to be deeper concepts, since they have straightforward generalizations to 
loop amplitudes [32], Amplitudes up to four loops have been constructed exhibiting the 
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duality and double copy [32, 40-51]. The implications of the color-kinematics duality have 
been used to derive a number of highly impressive results for string-theory amplitudes [71- 
78], and more generally the duality has been studied using string-theory methods [79-83]. 
Recently the double-copy construction has been extended to classical Kerr-Schild solutions 
in general relativity and gauge theory [84], 

In this paper the details of gravity amplitudes will not be discussed any further, since 
here we are instead interested in the aspect that the color-kinematics duality implies that 
primitive amplitudes obey the Bern-Carrasco-Johansson (BCJ) relations [31]. These rela¬ 
tions constrain the primitive amplitudes beyond the color-algebra basis discussed in the 
previous section. In the pure-gluon case these kinematic-algebra relations can be used to 
reduce the amplitude basis to (n — 3)! primitives. The simplest family of BCJ relations is 
linear in the generalized Mandelstam invariants: 4 

n —1 i 

X (X s -? n ) + 1,. • • ,n- 1) = 0, (2.15) 

i=2 j =2 

while the more general BCJ relations (found in ref. [31]) can be derived [85] from eq. (2.15) 
and its relabelings. 

The BCJ amplitude relations motivated the discovery of the scattering equations and 
associated elegant string-like formulas for gauge and gravity tree amplitudes [86-89]. In 
the context of scattering equations the BCJ relations have been applied to a massless 
quark-antiquark pair [53] and a massive pair of scalars [54], and three massive particles of 
different spins [55]. In section 4, we show that the color-kinematics duality (2.13) holds 
for general QCD amplitudes, including an arbitrary number of massive quarks, and give 
an explicit solution to the resulting BCJ amplitude relations. 

3 Color-algebra basis for quark-gluon amplitudes 

In this section we analyze the color structure of a general tree amplitude in QCD: we discuss 
the Melia basis [38] of (n—2)\/k\ color-ordered primitives and find their color coefficients. 
This will constitute a new color decomposition for QCD amplitudes. Though the results 
are applicable to arbitrary matter particles in any gauge-group representation, we refer to 
them as quarks. Similarly, the gluon can in principle be replaced by any adjoint particle 
since only group-theoretic properties are used in this section. 

We start with the case with no gluons, i.e. the multi-quark amplitudes, with all k 
quark lines having different flavors. In ref. [37] Melia considered such amplitudes with 
quarks in the adjoint representation and found a basis of primitives independent under all 
color-algebra amplitude relations. The latter are simply the familiar KK relations (2.10) 
projected onto the multi-quark case by setting to zero all the diagrams and primitives with 
crossed quark flavor lines (crossed with respect to a planar color ordering). 

This is possible because the KK relations can be deduced [31] from the fact that color- 
ordered primitive amplitudes can be expanded in terms of antisymmetric cubic vertices, as 

4 These relations for gluon amplitudes are sometimes called “fundamental BCJ relations” [85]. Here we 
avoid this terminology to prevent confusion with the fundamental representation of the gauge group. 
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is manifest in the adjoint representation, and can be imposed in the fundamental case (see 
figure 1). Under the projection that removes crossed fermion line diagrams, the antisym¬ 
metry of the vertices is maintained (since zero is antisymmetric) thus all the fermion-case 
KK relations are inherited from the pure-gluon case unscathed, though some of them are 
reduced to redundant equations, or even to trivial 0 = 0 equations. This implies that the 
fermion-case basis of primitives has to be a subset of the original basis of (n — 2)! primitives. 
In fact, the (n — 2)! basis decomposes into k\ = (n/2)! cases of inequivalent permutations 
of the quarks (maintaining the antiquark ordering), out of which only one has no flavor 
crossings. Hence, the basis surviving the projection is reduced in size by a factor l/(n/2)!. 
Indeed, the Melia basis for pure-fermion amplitudes has dimension (n — 2)!/(n/2)! [37], 
and in the general mixed quark-gluon case the dimension is (n — 2)\/k\ [38]. 

The fact that the basis of primitives were derived in refs. [37, 38] for the case of 
adjoint particles and color-ordered planar amplitudes is not a problem. The same planar 
amplitudes are also a basis of primitives in the mixed adjoint-fundamental case that we 
encounter in QCD. We illustrate this through examples, and then give the general color 
decomposition for QCD. 

3.1 Pure-quark example: n = 6, k = 3 

Let us now explain the details of the Melia basis using an instructive pure-quark example. 

The basis of primitive amplitudes for k = n/2 quark-antiquark pairs is given by [37] 5 

{-4(1,2, a) | cr € Dyck*,.! j , (3.1) 

where the last (n — 2) arguments must form a valid Dyck word out of the quarks and 
antiquark labels. 6 These words are defined as strings of X and Y letters, (k — 1) of each, 
such that the number of X’s preceding each Y is greater than the number of preceding 
Y’s. A more illuminating representation is obtained by realizing that Dyck words precisely 
correspond to well-formed brackets, with X’s playing the role of opening brackets and Y’s 
being closing brackets. For n = 6 there are two such words, XYXY and XXYY, equivalent 
to the brackets {}{} and {{}}, respectively. An example of an invalid Dyck word is XYYX, 
and it translates to the ill-formed brackets {}}{. 

Of course, the brackets themselves are not sufficient to specify the particle labels, in¬ 
stead they only specify the particle type. The X’s or opening brackets can be identified 
with quarks and Y’s or closing brackets with antiquarks (or, more precisely, the fundamen¬ 
tal and anti-fundamental representation, respectively). To form the Dyck words relevant 
for eq. (3.1), each valid bracket combination needs to be populated by particle labels. To 
be specific, if we have quark flavor lines 3-1—4 and 5 •(— 6, then there are exactly two label 
assignments per bracket combination that leave the flavor lines uncrossed: 

XYXY =>- (3,3, 5, 6 ), ( 5 , 6,3,4) O {3 4}{5 6}, {5 6}{3 4} , (3.2a) 

XXYY => (3,5, 6 , 4 ), ( 5 , 3, 4 , 6 ) O (3{5 6}4}, (5{3 4}6} . (3.2b) 

5 The precise basis used in ref. [37] is slightly different: {A(2,l, a) | a £ Dyck fe _ 1 } in our notation. 

6 Recall that we mark quarks and antiquarks with underscores and overscores, respectively. 
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5 ^ 3 . * P< g 5 < 13 ^ P< 


+ 



Figure 4: Feynman diagrams for the six-quark amplitude Ag r | e (l, 2,3,4, 5, 6). 


These four valid Dyck words are written out using two different notations: the bar notation 
used in eq. (2.11) and the bracket notation introduced above. As can be seen, the two 
notations convey the same information, since each downstairs and upstairs bar can be 
respectively identified with an opening and a closing bracket. In this paper we use both 
notations interchangeably according to convenience. A very important aspect of either 
notation is that each opening bracket corresponds to a unique closing bracket, and hence 
each downstairs bar corresponds to a unique upstairs bar. These pairs of states are precisely 
the quark-antiquark pairs of the same flavor. Thus the notation conveniently combines the 
information about the gauge-group representation and flavor. In other words, no separate 
notation is needed for specifying the flavor of the primitive amplitudes. 

We have thus arrived at the Melia basis for n = 6, k = 3, which contains four primitives: 

A(l,2,3,4,5,6), A(l, 2,5, 6 , 3,4), A( 1,2,3,5, 6 , 4) and A ( 1 , 2,5,3,4, 6 ) . (3.3) 

Now we will determine the color factors corresponding to the primitives and that are 
appropriate for quarks in the fundamental representation. Recall that the full amplitude 
can be written as 

4 tree _ c l n l , c 2 n 2 , C 3 n 3 ; C 4 U 4 ( ( CqTI§ ( C 77 I 7 ^ 

A ’ 3 - oT + o7 + oT + c7 + cy + cy + cT’ (3 ' 4) 

where each term corresponds to a cubic Feynman diagram in figure 4. Their color factors 


Cl 

_ rpa rpb rpa rpb 

~ 1 il3 1 j^2 1 hU 1 i^Q ’ 

_ rpb rpa rpa rpb 

c 2 — 1 i 1 j 1 ji2 1 izU 1 i^Q ’ 

(3.5a) 

C 3 

_ rpa rpb rpOu rpb 

~ 1 nj 1 j u 1 we 1 hi2 ’ 

r _ rpb rpa rpa rpb 

c 4 j i4, ± ibiQ ± il'i2 ’ 

(3.5b) 

C5 

_ rpa rpb rpa rpb 

~ 1 ib3 1 j^ 1 hi2 1 hu ’ 

r rpb rpa rpa rpb 

c 6 ± ie>j- L ji6 J -hi2 W4 ’ 

(3.5c) 


c 7 = -f abc 

rpa rpb rpc 

1 H12 1 hu 1 ’ 

(3.5d) 


can be read off according the rules in figure 1. For completeness and later use, we give the 
kinematic content of three representative diagrams: 

ni = -j(u 1 ^(l/i :3A +m 1 )Yv 2 )(u 3 'y ti V 4 )(u 5 'y u v 6 ), Di = (s lj3 , 4 - m 4 )s 34 s 56 , (3.6a) 

n 2 = ~j(uiY(¥i,5,6+rn 1 )'y ll V2)(u3'y t iV4 : )(u 5 'y I ,ve) , D 2 = (s 1)5 ,6 - rnl)s 34 ,s 5e , (3.6b) 

m = --(ui 7 m u 2 )(u 3 7 m u 4 )(u 5 (^i ) 2 - #3, a)v 6 ) + cyclic , D 7 = si 2 s 34 s 56 , (3.6c) 
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where we allow for arbitrary spins. All momenta are outgoing and the notation used is 


— ki + kj + • • • + ki , 

= (ki + % + ■■■ + &/) 2 , (3-7) 

Sij — {ki ~\~ kj) 

It is easy to see that the seven color factors (3.5) obey three commutation relations: 


Cl — C 2 = —C 7 , C 3 — C 4 = —C 7 , C 5 — C 6 = —C 7 . (3.8) 


Thus the color space of *4,0 r | e is four-dimensional, which is consistent with the length of the 
Melia basis (3.3). The primitive amplitudes therein can also be diagrammatically expanded 
according to the color-ordered Feynman rules (A.l) 


n 2 n 4 UR n 7 

l 123456 = JT 2 + ~Dl + 1 Tz + 1T 7 

_ ni , n 3 | n 5 n 7 

l I25634 - Dl + D 3 + D5 D 7 


™3 

77-4 

423564 “ jj 

~dI 

n 5 

n 6 

425346 — jj 

D 6 


(3.9) 


where for brevity we have written the amplitude arguments as subscripts. 

Now we wish to color-decompose the amplitude (3.4) such that we get exactly the Melia 
basis primitives as kinematic coefficients. We observe that we can massage it into a com¬ 
bination of the expressions in eq. (3.9) by using the commutation relations (3.8). Indeed, 
if we eliminate the color factors C 3 , cq and C 7 , we land on the following decomposition: 


4tree _ cini t c 2 n 2 , (ci - c 2 + c 4 )n 3 i c 4 n 4 | c 5 n 5 | (c 2 - ci + c 5 )n 6 , (c 2 - ci)n 7 

-~Eh + lh + W 3 + " d 7 + " d 7 + Di + ~D 7 

= C 2 ■'^423456 + Cl ^125634 + ( c 2 — c 4) ^123564 + (ci — C 5 ) ^^5346 ^ 

= 23456 ^123456 + ^425634 ^125634 + ^123564 ^-123564 + ^125346 ^125346 ■ 


The precise expressions for Cj^ ’s, the color coefficients of the primitives, are subject to 
the commutation identities (3.8). Our choice here is the one that will be generalized in 
section 3.3. Some traits of the general pattern can be seen from the following rendition of 
the color coefficients: 


C 


123456 


c 


125634 


■4 5 < 6 


5 123564 


■6 3- 


J ^125346 



(3.11) 


where we choose to draw some of the color diagrams in a nonplanar fashion in order to 
preserve the cyclic ordering dictated by the primitive amplitudes. In other words, although 
the color-ordered primitives (3.9) are composed only of planar Feynman diagrams, their 
color coefficients (3.11) contain both planar and nonplanar color diagrams. More than that, 
we observe that this non-planarity is related to the bracket “nestedness” for configurations 
{3{5 6}4} and |5{3 4}6}, as will be increasingly clear for higher-multiplicity examples. 
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+- 4 


A%f = 2 - 


3 + 4 


1 + 2 ' 


1+3 



5 5 

Figure 5: Feynman diagrams for the four-quark one-gluon amplitude -4!£| e (l, 2, 3,4, 5). 


3.2 Quark-gluon example: n = 5, k = 2 

Now let us replace the last quark pair by a single gluon, while four remaining quarks have 
the same distinct flavors as before. The quark masses are then mi = m 2 , m 3 = m 4 . For 
completeness and further use, we give full color and kinematic content for each of the five 


Feynman diagrams shown in hgure 5: 

ci = ThjTjiMsu > n i = 2^/2 ^ Ul ^ 1 , 5+m M tlv 2)(u3'y t iVi) , 

D\ = (S15-W?)S34 , 

(3.12a) 

„ _ rpb rjiCL5 rrib 

± ilj j 12 ± *3*4 5 

— ^ 

712 = 21/2 ( Ml7 ^ 2 > 5 rn ‘ 2 ^ 5V2 ')( u 3ln v 4), 

D 2 = (^25 1112)534 , 

(3.12b) 

_ pb pCL 5 pb 

6 ^*1*2 *3J ’ 

713 = 2^ ( Ml7 ^ V2 )(+^ 5 (^ 3 ’ 5+m3 ) 7/i,i;4 ) > 

D 3 = $ 12 (S35 m 2 ), 

(3.12c) 

_ rpb rpb rjiCl 4 

c 4 — - L iii2- L i33 ± ju 5 

—^ 

714 = 2 a /2 *- Mi7 ^ V2 ^ U37M ^ 4 ’ 5 TO 4)9 / 5f4), 

D 4 = S 12 (S45-ml) , 

(3.12d) 

„ _ ra^bcrj-ib pc 

C 5 — ./ ^2122^^3^4 ’ 

n 5 = ^((Ml95^2)(M3^5^4) - (U1^5«2)(W3^5V4) 

(3.12e) 


- (ui^ tt V 2 ){u 3 'y l j / V 4 : ){ki 2 -£ 5 )) > D 5 = S12S34 . 

The five color factors obey two commutation relations: 

ci — C 2 = —C 5 , C 3 — C 4 = C 5 . (3.13) 


implying a color-algebra basis of three primitive amplitudes. A symmetric color decompo¬ 
sition is obtained if we eliminate ci and C 4 : 


a tree _ 

5 ’ 2 ~ ^ Di 


= c 2 


n 1 


+ — 

Dl D 2 


n 3 n 4 \ 
D 3 Dj 


+ c 3 

= C 2 ^ 152.34 + C 3^-12354 + c 5^2 3 l! 




n4 U5 

+>4 ,0.5 


(3.14) 


12345 > 


where in the last line we identihed the correct color-ordered amplitudes, corresponding to 
the independent color factors C 2 , C 3 and C 5 (similar decompositions are found in ref. [90]). 
However, this is not yet the Melia basis [38]. To obtain it, we use the KK relation (2.10) 
to move leg 2 next to leg 1: 


A(1,J5 > 2, 3A ) = -Ml, 2,5,3,4) - A( 1 , 2,3,5,4) - A(l, 2,3,4,5) . (3.15) 

/3 OL 

Thus we obtain a new decomposition 

“^5,2 = —c 22li2534 — C 1 ^12345 + (~ c l + c 4)2li 


12354 


- Cl 2334 A-, nKJ >4 + C 1 n M t Hi n 315 + C 


(3.16) 


12534 ^-12534 ~ *- 12345 J H2345 ' '' 12354 ^12354 ) 
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in terms of three primitives 


A 


12534 


m 

d 2 


n3_ 

D 3 


. 

~A — “ . A 


n 5 

d 5 


_ m _ n 4 n 5 
M2345 - Di Da + D5 ’ 


A m* = w 3 + lk- (3I7) 


Note that they correspond to the single Dyck word XY, equivalent to the bracket {34}, 
with the gluon label 5 inserted before, after and in the middle of the word, respectively. 
Their color coefficients in the decomposition (3.16) are given by the following graphs: 


^12534 — 


^12354 — 




3 

2 




(3.18) 
4 , 

1 


where 6^2354 is drawn both as (— c\ + C 4 ) and (— c 2 + C 3 ) to emphasize that the nonplanar 
diagrams cannot be removed by commutation relations. The pattern to take note of is that 
the non-planarity occurs in the ordering {3 5 4} with the gluon sandwiched between the 
quark brackets, which is reminiscent of the nested quark-antiquark pairs in section 3.1. 


3.3 New color decomposition 

In this section we formulate the new color decomposition for QCD. 

We use the Melia basis for primitives with (n — 2k) gluons and k quark-lines [38]: 

{^4(1, 2, a) | cr £ Dyck fc _ 1 x {gluon insertions} n _ 2 fc} ■ (3.19) 

In the construction of this basis there are {2k—2)\/(k\(k — 1)!) Dyck words prior to assigning 
the particle labels inside the brackets. The quark labels can be assigned to {k — 1) slots in 
{k — 1)! inequivalent ways. The antiquark labels have unique slot assignments after this, 
since all quark lines have different ffavors. Then the (n — 2k) gluons are assigned to any 
place except between 1^ and 2, which must stay adjacent. With each gluon inserted, the 
number of available slots grows, starting from {2k — 1) up to (n — 2). Therefore, the size 
of this color-algebra basis is 

empty brackets 

x ' A * \ 

*{n,k)= x{k-l)\x {2k-l){2k)...{n-2) = ( -^^-, (3.20) 

v v y insertions of (n—2k) gluons 

dressed quark brackets 

in agreement with the reasoning given in the beginning of section 3. See table 2 for the 
explicit counts of the lower-multiplicity primitive amplitudes. 

Now the new color decomposition for QCD is conveniently written as 

x(n,fc) 

<?= E C{l,2,a)A{l,2,a), (3.21) 

crE Melia basis 
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k\n 

3 

4 

5 

6 

7 

8 

0 

1 

2 

6 

24 

120 

720 

1 

1 

2 

6 

24 

120 

720 

2 

- 

1 

3 

12 

60 

360 

3 

- 

- 

- 

4 

20 

120 

4 

- 

- 

- 

- 

- 

30 


Table 2: The number of n-point primitive amplitudes with k distinguishable quark pairs, inde¬ 
pendent under color-algebra relations, as given by the formula x(n, k) = (n — 2)\/k\. 


where A{\, 2, a) are the usual color-ordered planar primitive amplitudes (defined by the 
color-ordered Feynman rules in appendix A), and the color coefficients C(l, 2, a) are non¬ 
trivial objects that are the subject of the remaining discussion in this section. 

Using the suggestive bracket notation, we can obtain the color coefficients using the 
following replacement rules for the quark, antiquark and gluon labels: 


C(l,2,a) 


( — l) fc_1 {2|cr 11} 


q -> {,|T‘®S »_ 1 
9 ->• M 

n =“9 


(3.22) 


where the integer l is the level of bracket “nestedness” for a given particle in the word 
{2|ojl}. In other words, l is the number of opening brackets minus the number of closing 
brackets to the left of the particle. The bra {g| and the ket |g} now represent the funda¬ 
mental and anti-fundamental color indices of a quark and an antiquark, or can equivalently 
be understood as their the color wavefunctions. The object E? in eq. (3.22) is an operator 
obtained by tensoring l copies of the Lie algebra, 


"a _ 


El 


1(8>T a 


l® 1 


S=1 


(3.23) 


The sum effectively runs over each slot s in the identity tensor product and inserts a genera¬ 
tor with adjoint index a and in the appropriate representation, see figure 6. By convention, 
each copy of the Lie algebra representation corresponds to a particular nestedness level, 
starting from level l (the leftmost copy) and down to level one (the rightmost copy). Note 
that the operators S“ form a representation of the Lie algebra, 

[Sf.sf] =r lc 3f. (3.24) 


This makes it natural to extend eq. (3.23) to the pure-gluon case by defining Eg = T" d -. 

In the multi-sandwich formula (3.22) the quark and antiquark wavefunctions act only 
on the Lie algebra copy at their corresponding nestedness level l. For example, {21 and 11} 
act on the level-one copy of the group representation, which is complex conjugated with 
respect to the rest, as is seen in figure 6. This is indicated by the bar over the rightmost 
unit operator in eq. (3.23) and the fact that in eq. (3.22) we use the bra {2| and the ket 
11 } for the color wavefunctions of the antiquark 2 and quark 1, which is contrary to the 


15 - 




















Figure 6: Diagrammatic form of the operator E“. It is drawn as a single diagram with hollow 
quark-gluon vertices, this represents summation over the possible locations where the gluon line 
can attach. 


convention for all other particles. This is due to the special role of the fermion line 1-^—2. 
This notational subtlety could in principle be avoided by complex conjugating that line into 
1 —>• 2, i.e. by interchanging the roles of the quark and the antiquark.' More generally, any 
group representation copy can be complex conjugated, as it only changes the interpretation 
of which particles are quarks and antiquarks. (In fact, any representation can be allowed 
for each Lie algebra copy.) 

The above formulas (3.22) and (3.23) unambiguously determine the color coefficients 
in the color decomposition (3.21). For k = 1,0 it is straightforward to see that our de¬ 
composition coincides with the multi-peripheral formulas (2.11) and (2.12). Indeed, if all 
permuted particles 3, ...,n are gluons, they are replaced by E“* = T a \ producing pre¬ 
cisely eq. (2.11). In case particles 1 and 2 are gluons as well, the nestedness level is zero, 
and the adjoint-representation operators Eg 1 = yield the DDM decomposition (2.12). 

For a more detailed example of using eq. (3.22), let us consider the six-point amplitude 
discussed in section 3.1. One of the color coefficients in its decomposition is 

^123156 = {2|{3|T“® E?|4}{5|T fe ® 3?|6}|1} = {2|{3|T“® T a |4}{5|T b ® T b |6}|l} 

_. (3.25) 

= {2|T a r |l}{3|T a |4}{5|T fe |6} = {T b T a ) n - l2 T^ A T^ 6 , 

where we in the last step have translated it to a more conventional notation. A more 
interesting example is the color coefficient 

^23564 = (2|{3|T“<g) EU5|T b ® S||6}|4}|1} 

= {2|{3|T a ®T a {5|T fe (g)l®T ? '|6}|4}|l} + {2|{3|T a ® T a {5|T 6 ® T b ® 116}|4}|1} 

= {2|T a T b |l}{3|T a |4}{5|T b |6} + {2|T a |l}{3|T a T b |4}{5|T b |6} (3.26) 

_ / nnbrpa\ _ rpa rrib _ rpa (rparpb\ _ rpb 

1 1314-*■ Z 526 ^2112 v 1 i§lQ * 

The other two color coefficients for the six-point amplitude are, with less details, given by 
%5634 = {2|{5|T“®E((|6}{3|T 6 ®E?|4}|l} = {2|T a T b |l}{3|T f '|4}{5|T a |6}, 

Cmm = {2|{5| T a ® S?{3| T b ® S ft 2 |4}|6}|l} = {2\T a T b \l}{3\T b \4}{5\T a \6} (3.27) 

+ {2|T a |l}{3|T fc |4}{5|T a T 6 |6}. 

7 This would mean the Melia basis A(l,2, a), or, after relabeling, A(l,a,n). 
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The four above color coefficients are indeed the ones given diagrammatically in eq. (3.11). 

For the five-point amplitude discussed in section 3.2, our decomposition results in the 
following color coefficients: 

%53I = - {2|Hf{3|T fe <g> H$|4}|1} = - {2|T as {3|T b ® T b |4}|l} = - {2|T as T fc |l}{3|T 6 |4} , 

^12345 = ~{2\{3\T b ® Sj|4}S? 5 11} = -{2|{3|T b g>T b |4}T a5 |l} = - {2|T b T as |l}{3|T b |4} , 

%354 = -{2|{3|(T 6 ®S b 1 )S^|4}|l} (3.28) 

=-{2|{3|(T fe ®T 6 )(l(8)T a5 )|4}|l} - (2|{3|(T 6 ® T & )(T a5 <g) I)|4}|1} 

=-{2|T b T a5 |l}{3|T b |4} - {2|T b |l}{3|T b T a5 |4}, 

which coincide with the color diagrams in eq. (3.18). 

While we do not provide a proof for the decomposition (3.21), we have explicitly 
checked its validity for all quark-gluon configurations up to eight points, as well as the 
nine-point amplitude with four quark lines and one gluon. To be more precise, the check 
was done as follows. We expand the Melia basis of color-ordered primitives in kinematic 
cubic diagrams, 

A{1M= E ( 3 - 29 ) 

cr-color-ordered cubic graphs 

and solve for rii/Di in terms of the primitives A{\, 2, a). When this solution is plugged back 
into the color dressed cubic graph expansion (2.2), the color coefficients of the undetermined 
rii can be shown to vanish under the color algebra (2.4). What is left is then an expansion 
of the amplitude in terms of the primitives and some color factors c t . The color coefficients 
C(l, 2, <t) of the primitives can then be read off from this expression. They are indeed given 
exactly by eq. (3.22). 

3.4 Higher-point example and color-coefficient diagram 

Here we illustrate how the general formula (3.22) is applied to a high-multiplicity example. 
Consider the color factor of the following 14-point primitive amplitude with six quark pairs 
and two gluons (n = 14, k = 6): 

A(l, 2,13,3,5,6,4, 7,9,14,11,12, TO, 8). (3.30) 

We can use the cyclic property of planar amplitudes to move 1 to the end and then, to 
make the nestedness of the primitive more apparent, replace the bar notation of the legs 
with brackets: 

{2 13{3{5 6}4}{7{9 14{11 12}10}8}1} . (3.31) 

As before, the bar-bracket correspondence for the fermion line 1. <— 2 is opposite to all 
other legs due to the different convention for Lie algebra representation on that line. 

Using eq. (3.22) it is straightforward to obtain the expression for the color coefficient 
of this primitive 

^1,2,13,3,5,6,4,7,9,14,11,12,10,8 = — {2|.=. 1 13 {3|T fe (g) r, 1 {5|T (g) r, 2 |6}|4} (3.32) 

x{7|T d g)Sf{9|(T e (g)H|)S“ 14 {ll|T% -^112}110}18}11} . 
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Figure 7: Diagrammatic representation for the color coefficient of the planar amplitude 
A( 1,2,13, 3,5,6,4, 7,9.14,11,12,10.8), obtained by using the notation of figure 6. Note that the 
diagram has the same structure as the word {2 13{3{5 6}4}{7{9 14{1112}10}8}1}. 


While this formula is compact, it becomes a rather formidable expression if written 
in terms of standard color factors. In order to understand it better, let us sketch out the 
various standard color structures contained in it. For example, the first E“ 13 in eq. (3.32) 
will give rise to the structure 


{2|H“ 13 ...|l} = (T ai3 ...) l2 . i , (3.33) 

whereas the second appearance of S gives 

{2|...{3|T“®H?...|4}...|l} = (T“...). 3 _ 4 (...T a ...)_ 2 . i . (3.34) 


For l > 2 the operator Sf consists of a sum of tensor products, thus nested curly brackets 
imply summation over different possibilities of inserting the generator T a . For example, 
the third E in eq. (3.32) gives rise to two structures 


{2|{3... {5|2j® =||6}|4} ... |1} = (T% e (... (... 

+ (7*y. (...).- (...7* 

V '15*6 V '13U' 


) * 2*1 
)* 2*1 ’ 


while the sixth occurrence of E gives three contributions 


{2|... {7... {9 ... Eg 14 ... 110} 18}11} = (• ■ 

,. T aii ...)._ (...)._(...' 

/ 29U0 v ' *7^8 v 

^ *2*1 

+ (.. 

.). _ (... T“ 14 ). _(...' 

^ *2*1 

+ (•■ 

.). _ (•••).- (... T ai4 

^ *2*1 


(3.35) 


(3.36) 


Indeed, each S f gives rise to exactly l structures, which then are multiplied together. This 
implies that the number of standard color factors hiding in eq. (3.32) can be counted by 
multiplying the subscripts of the S’s in this expression, giving 


1x1x2x1x2x3x3 = 36 terms. (3.37) 

Finally, we note that probably the best way to understand the color coefficient (3.32) 
is to draw a diagram for it. Indeed, figure 7 contains the same information as the for¬ 
mula (3.32). In particular, compare this diagram with the word given in eq. (3.31). The 
diagram figure 7 is similar to a usual color factor diagram that describes the contractions 
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of f abc, s and T a ’s; however, we use hollow quark-gluon vertices to represent the operation 
of summing over the possible locations where a gluon can attach (see figure 6). In this sum, 
if the gluon does not attach at a given vertex, the vertex should be considered to be the 
identity operator. In this sense, the hollow vertex can either represent a gluon attachment, 
or a nonplanar crossing of the lines. Note that the diagram in figure 7, and similar ones 
built up by the diagrammatic H’s shown in figure 6, are the natural generalizations of the 
multi-peripheral diagrams that appear in the DDM decomposition. 

4 Kinematic-algebra basis for quark-gluon amplitudes 

In this section we turn our attention to the kinematic structure of general tree amplitudes 
in QCD. We find that they obey the color-kinematics duality (2.13), and this constrains 
the basis of primitives further. The number of independent primitives is reduced down to 
(n — 3)!(2 k — 2 )/k\ for k > 2, and for k = 1,0 it is given by the familiar (n — 3)!. 

The color-kinematics duality is well-established for gluons in any spacetime dimen¬ 
sion [31, 39, 46], and similarly for massless particles in supersymmetric YM multiplets 
that comprise gluons [32, 40-43, 45, 48, 51] or only matter [47, 49, 50]. The BCJ am¬ 
plitude relations, which are a consequence of the duality [31] but can be proven by other 
means [71, 72, 85], have in a few cases been shown to generalize to massive amplitudes. In 
particular, amplitudes with a massive scalar pair [54] and with three massive particles of 
different spins [55] have been considered in the literature. 

In this section, we show that the color-kinematics duality (2.13) holds for QCD tree- 
level n-point amplitudes with k quark-antiquark pairs having distinct flavors and masses. 
We prove this for k < 4, n < 8 by explicit calculations. Assuming the duality holds in 
general, we derive the resulting tree-amplitude relations. We expect that these relations 
not only hold for QCD, but can be applied to generic YM amplitudes with massive fla¬ 
vored matter particles (fermions/scalars), including general gauge-group representation, 
general spacetime dimension, with and without supersymmetry. In particular, the ampli¬ 
tude relations can be used as a gauge-invariant test to check if a given theory obeys the 
color-kinematics duality. 

We first give a few concrete examples and then state the general results. 
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4.1 Four-point example: n = 4, k = 1 


One of the simplest massive quark amplitudes involves one massive quark-antiquark pair 
and two gluons. It has three Feynman diagrams: 


3, a 4 ,b 



4, b 3,a 



4 ^2 (“193(^1,3+m)^2) = 

Z 513 — m 


rpb rpa 

7 /_//,/ \ / \ C2TI2 


2 S 14 —m 2 


Do 


(4.1a) 


(4.1b) 


3, a 


4 ,b 



■ fabcnnc 

13 (z{k A -Ez){ui¥4V2) - 2(k 3 -£ 4 ){uW3V2) 


2 S 12 


I,* 


2,1 


+(e 3 -£ 4 )(«i(#3 - 


(4.1c) 


£ 3^3 

D 3 


where the spins are left unspecified. These diagrams correspond precisely to the ones in 
the fundamental commutation relation (2.4b), which in terms of color factors read 


ci - c 2 = c 3 . (4.2) 

By the color-kinematics duality, the kinematic relation should then read 

ni — n 2 = n 3 . (4.3) 

Indeed, the kinematic numerators of the above Feynman diagrams, 

X X 

n\ = -- ui¥3(¥i,3+m)¥4V2 , n 2 = -- ui¥i(¥iA+ r n)fav 2 , 

. ^ ^ (4.4) 

n 3 = ^ 2 (fc 4 -£ 3 )(ixi^ 4 U 2 ) - 2 (fe 3 -e 4 )(tii 9 / 3 i; 2 ) + (e 3 -e 4 )(it 1 (^3 - # 4 )^ 2 )) , 

do satisfy eq. (4.3). To check this in detail one has to repeatedly use the Clifford algebra, 

t'V + 7 v = 27/^ , (4.5) 

the Dirac equations and gluon transversality conditions, 

ui(¥i — m) = 0, (^ 2 +m)u 2 = 0, ki ■ e, = 0, (4.6) 

as well as the mass-shell conditions for the quarks kf = /c 2 = m 2 and the gluons /c 2 = = 0. 

After several algebraic steps, one arrives at 

ni-n 2 -n 3 oc ui¥i ^¥^2 + u\ii¥$iV 2 - (e 3 ■ e 4 ) (^i (^1 + # 2 )^ 2 ) = 0 . (4.7) 

Note that the numerators (4.4) are gauge-dependent through the polarization vectors 
of the gluons, but the combination n\ — n 2 — n 3 is gauge-invariant (and zero). For a generic 
amplitude the numerators of cubic graphs need to absorb the four-gluon interactions, which 
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invariably leads to ambiguities in defining the numerators, and thus the duality is not 
manifest in general. However, in the above amplitude there is no quartic vertex, hence it 
is not surprising that the duality holds automatically from the Feynman rules. 

Having shown that the color-kinematics duality is present in this QCD amplitude, we 
proceed using only the formal properties of the numerators. Expressing the color-dressed 
amplitude through two color-ordered primitives is straightforward, 


A 


tree 

4,1 


3 


E 




ni n 3 
A + D 3 




C 2^1234 + Cl^4^243 > (4-8) 


where the color factors are of the multi-peripheral type, in accord with the decomposi¬ 
tion (2.11). Using the duality, n 3 = n\ — n 2 , we obtain a system of two equations, 


^1234 — 


/1 

1 > 

1 ni 

a ( 1 

1 ^ 

1 n 2 

[d 2 

+ Ac 

1 n2 D 3 

A I243 - [ D] 

+ Ay 

1 ni d 3 


which, after Gaussian elimination of m, yields 

A 


A- - ( J_ +J__ 

^1234 — l ^ w ^ 


Di +d 3 a ^ 3 ' 


(4.9) 


(4.10) 


,A D 3 ( A+A)A/ 

The coefficient of n 2 can be shown to be proportional to a vanishing sum of denominators 


A + A + A = (si 3 - m 2 ) + (si 4 - m 2 ) + s i2 = 0, (4.11) 


resulting in the following relation among the primitive amplitudes: 

(si 4 — 171 )j4i234 = ( S 13 — m ) A_243 ■ (4-12) 

This is a straightforward generalization 8 of the corresponding massless BCJ relation (2.15). 
It allows us to express the full amplitude in terms of a single primitive, 

= (igij- + A ma . <4.13) 

4.2 Kinematic algebra for n = 5, k = 2 

We proceed by showing that the color-kinematics duality is present in QCD amplitudes 
with two quark lines of different flavors and masses. Let us return to the n = 5, k = 2 
amplitude considered in section 3.2. Since the four-gluon vertex is also absent in this 
amplitude, the numerators in eq. (3.12) are directly given by the Feynman rules, and we 
expect the duality to hold automatically, 


ci — c 2 = — c 5 44 ni — n 2 = —n 5 , (4.14a) 

C3 — C4 = C5 44 U3 — U4 = TI5 . (4.14b) 

8 Eq. (4.12) can be presented [54] as (fci • ^ 4 ) 4^34 = (fci • ^ 3 )^ 4243 , making it formally identical to the 

massless BCJ relation, but here we prefer to use momentum invariants with explicit mass dependence, 

similar to Feynman propagators. 


- 21 - 






As before, to see that these relations indeed hold true, one has to use the Clifford algebra, 
the Dirac equation for massive quarks, the transversality condition for gluons, as well as 
the mass-shell conditions. We will spare the reader the details of this successful check, and 
instead proceed formally. 

Combining the two kinematic identities of eq. (4.14) and the three primitives as defined 
in eq. (3.14), we obtain a system of five equations. Four of these equations can be used to 
solve for numerators n\, n 2 , n 3 and in terms of variables 77 , 5 , A 15233 and A 12 354 , 


n 2 = ni - n 5 , 
n 3 = n 4 - n 5 , 


n 1 

714 


D1D2 

D\+D 2 

D3D4 

-D3+D4 


1^15234 ) > 


(A - ~ - —) 

1^12354 jj 3 > > 


the remaining fifth equation yields an expression for A 45345 , 


(4.15) 


D 2 Z?3 

A _—_ £ 4_ _ 4 —_u 

^12345 — n 1 n 15234 n i n ^12354 ^ 


D^D 2 ^ 


-D 3 + .D 4 


1 


D\+D 2 


+ 


1 1 

D 3 +D 4 


n 5 


(4.16) 


The variable 715 cannot be solved for since its prefactor in eq. (4.16) vanishes, resulting in 
the following gauge-invariant amplitude relation: 


(S 35 m 3)4^i2354 + ( s 12 S 34 )A ^ 2 345 ( s 25 m 2)^15234 — 0 . (4-17) 

This expression involves the primitive A 152 34 outside of the Melia basis, but we can trade 
it for A 1253 4 using the KK relation (3.15). Via the identity S 12 — S 34 = S 35 + S 45 — m| — m|, 
we obtain another amplitude relation: 

(s 25 ~ rn 2 2 )A^ 34 + (s 2 5 + S 35 — — m3)A i2 3 53 1 

+ (S 25 + S 35 + S 45 m\)A-fi^ = 0. 

Here, in order to have consistent quark masses, one should set mi = m 2 and m 3 = 1714 . 
We keep the masses distinct, as well as avoid using further momentum identities, so as to 
match the pattern of the general n-point formula given in the next section. If we do use 
further momentum identities this expression can be simplified to 


( s 25 m 2)^-12534 4“ ( Sl4 s 23)A3 2 3 53 (^15 m l)^l_2345 — (4.19) 


which is nothing but a relabeling of eq. (4.17): (1,2,3,4) —> (4,1,2,3). In fact, all BCJ 
relations for n = 5, k = 2 are related by such simple relabelings. 

As an alternative test of the color-kinematics duality, one can verify these amplitude 
relations directly, thus circumventing the numerator check. For example, it is straightfor¬ 
ward to check that the relations (4.17) and (4.19) hold in the m; —» 0 limit, in which the 
amplitudes become 9 


A(r,5 + ,2 + ,3-,4 + ) 

A(r,2 + ,3-,5 + ,4 + ) 


■ < 13 > 2 

(15) (25) (34) ’ 

■ < 13 > 2 

(12) (35) (45) ’ 


A(r,2 + ,5 + ,3-,4 + ) 

A(r,2 + ,3-,4 + ,5+) 


■ (13) 2 (23) 

(12) (34) (25) (35) ’ 

. (13) 2 (14) 

(12)(34)(15)(45) ’ 


9 Since k > 1 these primitives do not correspond to component amplitudes in any of the pure Af = 1, 2, 4 

SYM theories. Thus their BCJ relations in principle need not be the same as those of pure SYM theories. 
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where we have specialized to quarks and a gluon of definite helicity, using the spinor- 
helicity formalism [91]. Of course, the above amplitude relations hold just as well for any 
mass and spin/helicity configuration, as follows from the fact that the Feynman diagrams 
in eq. (3.12) satisfy the color-kinematics duality. 

We note that eq. (4.18) can be mapped to the gluonic BCJ relation (2.15) by letting 
Sj 5 —>■ ( s j 5 ~ m]) = 2 kj ■ k§. However, not all of the gluonic BCJ relations hold true for 
quark primitives. For example, consider the following identity valid for gluons [31]: 


114235 — 


s 12 545^412345 + S25(si4 + s 24)^-15234 


(4.21) 


S 35 S 24 

In trying to map the primitives to the quark flavor/mass configurations considered here, 
the left-hand side involves crossed flavor lines and thus should be set to zero. The resulting 
two-term relation between ^5345 and A 15 2 3 4 is not a true identity, not even in the massless 
limit. This shows that mapping the gluonic BCJ relations onto the k > 2 quark case is 
more subtle than for the corresponding KK relations, which we argued in section 3 to hold 
true regardless of the particle configuration. Nevertheless, it is not difficult to understand 
the general structure of the quark-gluon BCJ relations, as we show in section 4.4. 

Before proceeding to higher-point amplitudes, we note that one can use eq. (4.18) to 
express the full amplitude (3.16) in terms of two primitives with legs 1, 2 and 3 fixed: 


4 tree _ 
^5,2 — 


rjib rpCL 5 rpb . rjib nnOb§ rpi 

1 iii2 1 iz3 1 ju ' 1 hj 1 ji2 1 hu 


b S 35 - m 3 


S 25 - 


( rjids rpb rpb 1 rpb rpCl^ rpb '^15 TTl^ 
1 il3 1 ji2 1 hu ' 1 hj 1 ji2 1 i3U s ^ _ m 2 


^12354 


^12345 


(4.22) 


4.3 Kinematic algebra for n = 6, k = 3 


Let us now consider the kinematic structure of the six-quark example of section 3.1. The 
general massive forms for three of the numerators are found in eq. (3.6), and in the limit 
of massless quarks these simplify to 


ni(r,2 + ,3-,4 + ,5-,6 + ) = —i<13)(5|l+3|4][62] , 
n 2 (l“ 2 + ,3-,4 + ,5",6 + ) = -i(15)(3|l+5|6][42] , 

77 , 7 (1 ”, 2 + , 3“,4 + , 5“, 6 + ) = —z(<13)[42]<5| 1+2|6] + <35>[64](1|3+4|2] 

+ (15)[62](3|5+6|4]) . 


(4.23a) 

(4.23b) 

(4.23c) 


The remaining four numerators are obtained by permuting the particle labels of these 
expressions. As is evident from this form, it so happens that all the Feynman diagrams are 
separately gauge-invariant. Indeed, not only is the quartic vertex absent, but furthermore 
there are no external gluons to introduce gauge dependence (e.g. via a reference axial 
vector), as it happens in the five-point amplitude. Therefore the numerators are unique, 
and hence the color-kinematics duality cannot be imposed unless it holds from the start: 


Cl - C 2 = c 7 , 
C 3 — C 4 = —C 7 , 
C 5 — C 6 = —C 7 , 


n\ — n 2 = —n 7 , 
n 3 — n 4 = —n 7 , 
n 5 -n 6 = -n 7 . 


(4.24) 
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Fortunately, these numerator identities do hold, as can be checked for the massive numera¬ 
tors (3.6) through Dirac spinor algebra, or in the massless limit (4.23) using spinor-helicity 
identities. 

If we now return to the formal definitions of the primitives (3.9), and consider them 
together with the kinematic algebra of eq. (4.24), we obtain a system of seven equations 
for seven numerator variables, which is straightforward to solve uniquely. The first six 
numerators are given by 


n i = 


M( 


l 


+ 


l 


-D 3 +D 4 D 5 +D e d 


j ) n 7 ^125634 + 


d a 


-D3 + D4 


-A 


123564 


+ 


D 6 


D^+Dq 


-A 


125346 


}> 


ri 2 =ni+n 7 , 
n 4 = n 3 + n 7 , 


Ds 

ns d 3 +d 4 
d 5 

125 D 5 +D 6 


(^ 4^123564 

( Dq A . jj 25346 


+ n 7 ), 
+ n 7 ), 


n 6 = n 5 + n 7 , 


(4.25) 


and the seventh numerator is 


n 7 


|(D3 + -D4)(D5 + Z?6)(-D2j4jj23456 -^1^125634) {kD 3 +Dq)(D 2 D 3 D 1 D 4 ) 

~(D 3 + D 4 )(D 2 D 5 - D[ 

/{{D 1 +D 2 )(D 3 + D A ){D 5 +D 6 )/D r - (D 3 +D A )(D 5 +D 6 ) 

— (Di+D 2 )(D 3 + D4+D 3 + Dq)^ . 


(4.26) 


Indeed, we have now expressed all the m in terms of gauge-invariant quantities, formally 
confirming the direct observation that each Feynman diagram is gauge-invariant. The 
above expressions look non-local, but by uniqueness of the solution, the numerators have 
to be local, as is manifest in eqs. (3.6) and (4.23). 

An interesting fact we learn from this solution is that there are no extra relations gen¬ 
erated 10 among the primitive amplitudes in the color-algebra basis (3.9). More generally, 
this curious absence of BCJ relations happens when there are no external gluons, n = 2k, 
as we have explicitly checked up to eight points. At eight points this fact is somewhat 
counterintuitive since the Feynman graph expansion involves the quartic gluon vertex that 
naively should lead to numerator ambiguities and associated BCJ relations. 

4.4 BCJ relations for QCD 

In the previous sections we exemplified that scattering amplitudes in QCD obey the color- 
kinematics duality. For these examples, with many quarks and no quartic gluon vertices, 
it turned out that the color-kinematics duality followed directly from the Feynman rules. 
However, in general this is not the case, and the duality has to be imposed by hand 
using nontrivial rearrangements of terms between the kinematic numerators. This is what 
happens for higher-point amplitudes with only gluons [31] or with one quark-antiquark 
pair (n > 5 ,k < 1). For such amplutudes the numerators satisfying the color-kinematics 

10 This does not imply that further amplitude relations do not exist. However, one should expect that any 
extra relations require more that just the color properties (2.3), (2.4) and their kinematic analogue (2.13). 
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duality were shown to exist to all multiplicities [80, 92, 93], thus effectively proving the 
duality at tree level for k < 1. 

By explicit calculations we have checked that the duality works for any quark-gluon 
configurations up to eight particles. While we do not provide a proof for n > 8, k > 1, 
we will show that the duality imposes the BCJ relations in QCD that are constitute a 
well-defined subset of the pure-gluon BCJ relations. 

In the pure-gluon case, the number of independent BCJ relations are (n — 2)! — (n — 3)! 
and they are given in ref. [31]. A subset of those are linear in momentum invariants, 

n— 1 i 

^2 + = 0, (4.27) 

and in ref. [85] it was shown that relabelings of these simple equations could be used to 
derive the more complicated relations involving higher powers of momentum invariants. 

We found that the corresponding quark-gluon BCJ relations, for k massive quark lines, 
are given by the general formula 

n —1 i 

y; (y s in — m 2 }y4(l, 2,... i, to, i + 1,... ,n — 1) = 0 . (4.28) 

i= 2 j= 2 

where particle to is strictly a gluon, while the remaining (to — 1) particles can be of any 
type: quark/antiquark/gluon. In the next section we will generalize this formula to include 
the relations with higher powers of momentum invariants, as well as derive the number of 
linearly independent BCJ relations, counted for n < 8 in table 3. 

In section 4.2 we derived a to = 5, k = 2 amplitude relation (4.18) with the permuted 
leg to being a gluon. It has precisely the form (4.28). The four-point relation (4.12) from 
section 4.1 does not have this precise form right away, it can be easily rewritten that way, 
either as a sum over different insertions of the gluon leg 4, 

(S24 ~ m ^)^1243 + ( S 24 + S34 — m2 ) A l234 = 0 , (4.29) 


or, equivalently, as a sum over insertions of the gluon leg 3, 


(s 2 3 - rTO 2 )A i 234 + (s 2 3 + S 3 4 ” m2 ) A 1243 = 0 • 


(4.30) 


For the pure-quark six-point amplitude in section 4.3, we note that eq. (4.28) is consis¬ 
tent with the fact it had no BCJ amplitude relations despite obeying the color-kinematics 
duality. Indeed, by definition pure-quark amplitudes, to = 2k, have no external gluons, and 
thus eq. (4.28) gives no relations for them. 

We derived the quark-gluon BCJ relations (4.28) as follows. We start with the primitive 
amplitudes in the Melia basis (3.19), 


{-4(1,2, a) y ± 

cr-color-ordered cubic graphs Tj 


a £ Melia basis 


is} , 


(4.31) 
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5(n, k) 


k\n 

3 

4 

5 

6 

7 

8 

0 

0 

1 

4 

18 

96 

600 

1 

0 

1 

4 

18 

96 

600 

2 

- 

0 

1 

6 

36 

240 

3 

- 

- 

- 

0 

4 

40 

4 

- 

- 

- 

- 

- 

0 


(n — 2)! — (n — 3)! 

((n- 2)! - (n- 3)!(2fe - 2 ))/k\ 


for k = 0,1 
for 2 < 2 k < n 


Table 3: Number of independent BCJ relations, S(n, k), for n-point amplitudes with k distin¬ 
guishable quark pairs and (n — 2k) gluons. 


written in terms of the formal numerators rtj, and the given kinematic denominators D j. 
The numerators are further constrained by the kinematic Jacobi/commutation relations: 11 


rii - rij = n k 


(i,j,k) € Jacobi/commutation triplets j 


(4.32) 


Solving the combination of the systems (4.31) and (4.32), similar to how it was done in 
sections 4.1-4.3, one obtains a collection of gauge-dependent solutions for a subset of the n* 
and a collection of gauge-invariant solutions for a subset of primitives. The latter solutions 
are the BCJ relations. 

We explicitly solved this system for any quark-gluon configuration up to eight particles, 
as well as for the nine-point amplitude with four quark lines and one gluon. We also 
verified that, similarly to the pure-gluon case, all the BCJ relations in QCD follow from 
label permutations of the simple family of relations of the type (4.28). However, the label 
permutations of eq. (4.28) should be regarded as an equation system, not a solution, to 
the full set of BCJ relations. In the full solution the primitive amplitudes are reduced to 
a new basis smaller than that of the color-algebra basis (3.19). 


4.5 New amplitude basis for QCD 

Now let us find a complete solution to all gauge-invariant relations imposed by the system 
of equations comprised of eqs. (4.31) and (4.32). These relations are the BCJ relations and 
allow us to reduce the primitive amplutudes to a smaller basis. 

The new BCJ basis for k > 2 is taken to be a subset of the Melia basis (3.19) obtained 
by demanding that the third particle is strictly a quark: 

{HQ, 2, g, a) | { q,a} € Dyck fc-1 x {gluon insertions in cr} n _ 2 fc} . (4.33) 

It is always possible to demand that for k > 1, because any gluon that happens to be in 
the third location can be moved to other positions using eq. (4.28). 

Let us count the size of the basis (4.33). As before, there are (2k — 2 )\/(k\(k — 1)!) 
Dyck words, and after dressing them with quark labels, each allow for (k — 1)! different 

n The number of independent Jacobi/commutation relations is u(n, k) — >c(n, k) = — Lirpi, 



















P(n,k) 


k\n 

3 

4 

5 

6 

7 

8 

0 

1 

1 

2 

6 

24 

120 

1 

1 

1 

2 

6 

24 

120 

2 

- 

1 

2 

6 

24 

120 
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- 

- 

- 

4 

16 

80 

4 

- 

- 

- 

- 

- 

30 


(n — 3)! for k = 0,1 

(n — 3)\(2k — 2)/k\ for 2 < 2k < n 

Table 4: Number of independent primitive amplitudes, /3(n, fc), in the full n-point amplitude with 
k distinguishable quark pairs and (n — 2k) gluons, after imposing the BCJ relations. 

quark pair configurations. The (n — 2k) gluons are then free to be assigned to the slots 
in-between the particles, except for the space inside the fixed sequence 1, 2, q. For each 
gluon inserted the available number of slots increases, ranging from (2 k — 2) to (n — 3). 
Thus the length of the basis is 

empty brackets 

/ A \ 

P(n,k)= ^ _ 2 ^’ ! x(fc - 1)! x (2fc - 2)(2 k - 1)... (re - 3\ 

^ ' v ** insertions of (n—2k) gluons 

dressed quark brackets 

For the amplitudes that only have one quark/antiquark pair, k = 1. it is not possible 
to pick the third particle to be a quark. So instead we pick it to be the gluon 3, giving 
the basis A(l, 2, 3, <r). Except for the bars on the labels, this is the same basis as for the 
pure-gluon case [31], thus the basis for k = 0,1 is of size (n — 3)!. The basis counts for 
different quark-gluon configurations are exemplified and summarized in table 4. 

The full solution to the BCJ relations are given by first solving the numerators rq in 
terms of the primitives in the BCJ basis (4.33), and then plugging them into the primitives 
that are not part of this basis. Since the color-algebra basis (3.19) already has legs 1 and 2 
next to each other, we only need to give the reduction formula for the primitives with the 
first quark leg q separated from leg 2 by a set of gluonic legs a and followed by a mixed 
set of quark-gluon legs (5: 

A(l,?,a,q,P). (4.35) 

To simplify the subsequent formulas, we choose q = 3 and take the leg labels in sets a and 
/3 to be consecutive numbers: 


(n — 3)l(2k - 2) 
k\ 


(4.34) 



a = {4,5,... ,p — l,p} , q = 3, /3 = {p + l,p + 2, ..., n - 1, n} . (4.36) 

As already mentioned, a consists strictly of gluon legs, and the particles in f3 can be of 
any type: quark/antiquark/gluon. The consecutive labeling choice can always be undone 
in the final expressions by relabeling of legs 3,4 ,,n. 




















By extrapolating from the structure of the solutions evaluated up through eight exter¬ 
nal particles, we obtain the following all-multiplicity formula: 


A(l, 2 , a, 3, j3) 


= Y, 2, 3, a) 


a-eS'(a)LU/3 i =4 


CJ, 1 | i) 

S2,4,...,i - ’ 


(4.37) 


where the sum runs over partially ordered permutations of the merged sets a. and (3. The 
shuffle product is used in S(a ) ULI/3 to obtain all permutations of a U (3 that preserve 
the order of the elements of (3 . Note that such permutations maintain the valid bracket 
structure of the quark labels. The kinematic function T associated with the gluon leg i is 
given by 


J r (3,(Ti,(T 2l .. • ,o- n _ 3 , l\i) = F(p\i) 

(4.38) 


and where tk is the position of leg k in the set p, except for 1 3 and t p+ 1 , which are always 
defined to be 12 

t'i = ^5 , tp +1 = 0. (4.39) 

(Note that for p = 4 this implies 1 3 = t p+ 1 = 0.) The function Q is given by 

re ) s v ~ m i if * < 3 or 3 = 1 : 3 1 . . 

5 («) = | 0 else )• (« 0 ) 


if ti-l < ti 
-YdLi 0{i,Pi) if ti-i > ti 


+ 


®2,4,...,i - ml if <ti < t i+ 1 

®2,4,...,j + ml if U- 1 > tj > tj+i 

0 else 


For the case that a consists of a single gluon, the formula (4.37) becomes equivalent to 
eq. (4.28). If the amplitude has a single massive quark line, k = 1, the relations (4.37) still 
hold after converting leg q to a gluon: 3 —>■ 3. Similarly, in the pure-gluon case, dropping 
all bars reduces eq. (4.37) to the standard BCJ relations [31]. We have explicitly checked 
that all QCD amplitudes through eight points satisfy the BCJ relations (4.37). 

A known feature of the BCJ relations is that they are very general. We expect eq. (4.37) 
to be valid not only for ordinary QCD, but for more generic gauge theories describing the 
interactions of adjoint vectors and massive matter transforming in any representation of 
the gauge group. In particular, supersymmetric versions of QCD, as well as D-dimensional 
extensions of QCD, are examples of theories that should obey eq. (4.37). At tree level there 
is little difference between super-QCD amplitudes and those of QCD. The reason is that the 
gluons and quarks directly translate to supersymmetric adjoint and fundamental multiplets, 
and the structure of the supersymmetric tree-level interactions are dictated by those of the 
gluon and quark states (which are top components of respective supermultiplets). 

The pure-gluon BCJ relations are known to be valid in any space dimension [31, 71, 
72, 85], and we find the same is true for the quark-gluon amplitudes considered in this 
paper. For example, we used only D-dimensional properties to show that the numerators 

12 An alternative choice is ts = oo,t p +i = 0, which is equivalent to eq. (4.39) by momentum conservation. 
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generated from the Feynman rules in sections 4.1-4.3 obey the color-kinematics duality. 
The observation that it holds for D-dimensional quark-gluon amplitudes explains why the 
generalization to massive (Dirac) quarks is straightforward. As is well known, the mass of 
the quarks can be reinterpreted as a component of a higher-dimensional massless momenta, 
living in some an internal space (i.e. similar to the Kaluza-Klein construction). Indeed, the 
factors (s 2 , 4 — m|) and — mj) that appear in eqns. (4.37-4.40) can be thought of as 
being higher-dimensional versions of S 2 , 4 ,...,j and S{j , given that the quarks have momenta 
antiquarks momenta (fcf,— Wj), and gluons have momenta (fcf, 0). 


4.6 New mixed decomposition 

Finally, we can combine our color decomposition (2.9) with the solution to the BCJ relations 
to decompose the full amplitude in the BCJ basis of primitives (4.33). 

For simplicity, let us start with the case of k = 0,1. We can rewrite both Eqs. (2.11) 
and (2.12) as 

<| e = E C(l,2,a)A(l,2,a), (4.41) 

<r£Sn-2({3,...,n}) 

where <7(1, 2, a) = (T a -W ... T 0 ^")) in the fundamental or adjoint representation. The 
sum over permutations of arguments 3,..., (n — 3) can be divided into those terms with 
with particles 1, 2, 3 appear together, and those terms where 2 and 3 are separated by a, 


<?= E <7(1,2,3,<4.4(1,2,3,cr) + £ <7(1, 2, a, 3, /3) A(l, 2, a, 3, f5 ). (4.42) 

o-€S„_3({4,...,n}) «U/3={4,...,n) 

an/3= 0 

The BCJ relations [31] allow us to convert all the primitives in the second sum to the BCJ 
basis, with the three first legs fixed, 


A(l,2,a,3,/3) 


= E 2 ’ 3 , *) n 

crG5'(a)LLI/5 i —1 


F(3, cr, l|i) 

■ 52 , 01 , ...,Q!i ^2 


(4.43) 


where for non-canonical orderings the function T is defined via label permutations of 
eq. (4.38). By construction, since the permutation o is an element in S(a) LU/3, it follows 
that (3 is a subset of a. Hence, after plugging eq. (4.43) into eq. (4.42), one can rearrange the 
order of summations by defining f3 C cr and a € S(cr\/3). This gives the BCJ decomposition 


A 


tree 

n,k< 1 


E ( 4 - 44 ) 

oG5 ri _3({4,...,n}) 


X 



3,011 

2—1 


F{2>,crA\i) \ 

S2, ai ,..., ai - ml) ' 


This expresses the full amplitude in terms of the basis of primitives independent under 
both KK and BCJ relations. Note that their gauge-invariant coefficients now contain not 
only color but also kinematic factors. 

It is now straightforward to generalize the decomposition (4.44) to the general quark- 
gluon case. For that, the external sum should run over the BCJ basis (4.33), in which the 


- 29 - 




quarks q directly succeed 1, 2, while the internal sum over (3 should include all quark and 
antiquark labels, and thus the complement set cr \ f3 should be purely gluonic: 


<?>2= E 


(4.45) 


(<?,cr) EBCJ basis 


M 


C(l,2,q,a) + E E C(l,2,a,q,P)l 


l\i) 


(3Ccr aeS(cr\(3) 
a\/3 gluonic 


^ s 2,ai,...,ai 


where the color structures are defined as in eq. (3.22). This decomposition should hold for 
any k = 2,..., \n/2\. In fact, it also holds in the k = 0,1 case. Indeed, once the bars are 
dropped and q —>• 3, it becomes indistinguishable from eq. (4.44). In the pure-quark case, 
n = 2k, the internal sum over the gluon set a vanishes, and the BCJ basis collapses to the 
Melia basis, making eq. (4.45) equivalent to the color decomposition (3.21). 

We have explicitly checked the decompositions (4.44) and (4.45) through eight points, 
as well as for the nine-point case with four quark-antiquark pairs. 


5 Summary and discussion 

In this paper we explored and organized the color and kinematic content of general tree 
amplitudes in QCD with flavored massive quarks and massless gluons. We decomposed 
the amplitudes into reduced sets of color-ordered primitive amplitudes and associated color 
coefficients, and found that the amplitudes obey the color-kinematics duality. 

A familiar observation is that the true color space of a gauge-theory amplitude is 
smaller than the space spanned by the classic SU(IV C ) color basis expansion. Exploiting 
this fact, we obtained a new color decomposition for n-point QCD amplitudes with k quark- 
antiquark pairs, which involves only the (n — 2)\/k\ planar color-stripped primitives that 
belong to the Melia basis [37, 38]. This decomposition can be regarded as a natural (n, k)- 
generalization of the known Del Duca-Dixon-Maltoni decomposition for gluons [33, 34], 
which uses only the (n —2)! color-ordered primitives comprising the Kleiss-Kuijf basis [30]. 

The DDM decomposition is advantageous because of its analytic compactness. For ex¬ 
ample, it has been successfully used for efficient computations of the color-subleading parts 
of QCD loop amplitudes [65, 94], Similarly, one can hope that the new decomposition 
can be used to improve on the efficiency of calculations relevant for LHC phenomenol¬ 
ogy. However, for calculations of color-averaged cross-sections, properties other than those 
sought after in formal amplitude calculations may become more important. For instance, 
it can be more convenient to have color factors that are orthogonal to each other (e.g. see 
the recent ref. [95]). The color factors used in the decomposition presented here consist 
of those appearing in the Feynman diagrams of the amplitude, thus they are completely 
general and can encode any gauge group and any group representation. 

The new color decomposition has an interesting hierarchy of nested structures of quark 
lines. For the least nested configurations the color factors are planar and similar to the 
multi-peripheral diagrams, whereas the color coefficients of the more nested primitives have 
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an intricate nonplanar (with respect to particle ordering) structure. This structure is con¬ 
trolled by multi-tensor representations of the gauge group, with generators suggesting 
that a deeper understanding of the Lie-algebraic structure of the decomposition could be 
advantageous. Using the higher-representation generators S“, we obtain surprisingly com¬ 
pact expressions for the color coefficients of the primitive amplitudes, which can be rather 
formidable expressions when expressed in terms of ordinary color factors. The color dia¬ 
grams that we use for the color coefficients generalize the multi-peripheral diagrams used 
in the DDM decomposition. 

Nested color structures of quark lines and gluons also appear in other gauge-invariant 
observables, e.g. in processes where massive quarks and hard gluons can be approximated 
by Wilson lines. These are important for understanding the soft and collinear singularity 
structure of gauge-theory amplitudes [96]. Non-abelian exponentiation theorems allow the 
Wilson line calculations to be reduced to simpler web diagrams. Similar to the amplitudes 
considered in this paper, the color factors arising from webs are not independent, and the 
mixing of color and kinematic degrees of freedom, which plays an essential role for the 
infrared exponentiation [97-99], is governed by nontrivial combinatorics [100]. It would be 
interesting to better understand the common aspects of the new color decomposition and 
well-known properties of Wilson lines and webs. 

In this paper we provided substantial evidence that the color-kinematics duality [31, 
32] is present in QCD. In particular, we considered the duality for the first time in the 
context of arbitrary numbers of distinctly flavored massive fundamental Dirac fermions 
(quarks), and show that they mesh well with it. By explicit calculations up to eight 
particles (quarks and gluons), we showed that one can find kinematic numerators that obey 
the same commutation/Jacobi identities as the corresponding color factors. This suggests 
that a kinematic algebra, analogous to the color algebra, controls the QCD amplitudes. 
We do not explicitly consider loop amplitudes in this paper, but given the existence of 
the duality at tree level we conjecture that it should also be present in loop amplitudes 
of QCD. Indeed, through the unitarity method [3-5] we know that any unitarity cut, 
which breaks up a loop amplitude into products of tree amplitudes, will inherit the color- 
kinematics duality of the corresponding trees. Every such unitarity cut should then have a 
diagrammatic expansion where the numerators obey the color-kinematics duality, which is 
highly suggestive of a kinematic algebra that exists at loop level. Certain simple massless 
QCD amplitudes have already been shown to obey the color-kinematics duality at one and 
two loops [44-47, 50], adding to the credibility of the conjecture. 

Using the color-kinematics duality, we derived a complete set of gauge-invariant BCJ 
amplitude relations for QCD at tree level. These amplitude relations are massive (n, k)- 
generalizations of the purely gluonic BCJ relations [31]. While the analytic form of the 
new relations is similar to the pure-gluon relations, their implications differ substantially. 
In particular, the new relations clarify the different roles that gluons (vectors) and quarks 
(matter) play in the context of the color-kinematics duality. The existence of the BCJ 
relations is directly tied to the presence of external gluons. Pure-quark amplitudes have no 
BCJ relations, while pure-gluon amplitudes have the maximum number of such relations. 
The most elementary type of BCJ relations corresponds to the action of moving a single 
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gluon around and inserting it at possible locations along the color-ordered amplitude. This 
is similar to the action of a Ward identity, or, indeed, the action of a Lie algebra. That 
BCJ relations arise only from the presence of external gluons supports the observation that 
these relations are intimately tied to the gauge freedom of the theory [31, 32, 39]. 

Using the new BCJ amplitude relations, we found a new basis of (n — 3)!(2 k — 2 )/k\ 
primitives for k > 2, which should be contrasted to the well-known BCJ basis of size 
(n — 3)! for pure-gluon or single-quark-line amplitudes, k = 0,1. The new basis allows 
us to construct an amplitude decomposition that uses only these primitives. In doing so, 
we introduce a mixture of kinematic dependence and dependence on color factors into the 
coefficients of the primitives. It is interesting to note that for k > 2 this decomposition 
sidesteps (n — 2 k) (n — 3) \/k\ of the Melia basis primitives, and for k = 0, 1 it sidesteps (n — 
3)(n —3)! of the Kleiss-Kuijf basis primitives. For the reader’s convenience, a Mathematica 
implementation of the amplitude decompositions is provided in the ancillary hie [101]. 

Finally, we note that due to the duality of color and kinematics one should be able to 
swap color structures for kinematic structures in several formulas of this paper, since they 
obey the same general Lie-algebraic relations, from which most results were derived. For 
example, if one takes the new color decomposition introduced in section 3.3 and replaces 
the color coefficients C(1,2 ,cj) with the corresponding kinematic coefficients K(l,2,a), it 
should give a gravitational scattering amplitude 

M%£= £ K(l,2,a)A(l,2,a), (5.1) 

crEMelia basis 

where the gravitational coupling has been suppressed. Analogous to the color coefficient, 
K(l,2,cr) is a local kinematic function that consists of a sum over kinematic numerators. 
Like the numerators, this function is generically gauge-dependent. By the color-kinematics 
duality, K (1, 2, a) should have a formula that mimics the formula for C(l,2,a). This sug¬ 
gest that the higher-representation generator Sf should have a kinematic analogue, which 
would represent the action of a gluon on tensor products of fermion lines. Improved under¬ 
standing of this kinematic object would most likely be a concrete step towards unraveling 
the kinematic Lie algebra of QCD. 
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A Color-ordered Feynman rules 


For completeness, in this appendix we give the color-stripped Feynman rules [56], which 
are consistent with the color vertices in figure 1. All momenta are considered outgoing. 
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